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ABSTRACT 



Satel 1 i te 


detection 


of 


complex aerosols. 


in 


particular, oil 


smoke over 


water 


and dust over land. 


is 



generally difficult. On 1 March 1991, a smoke plume generated 
by burning Kuwaiti oil wells and a dust storm over the 
southwestern Arabian Peninsula, provide the opportunity to 
study both of these effects. Utilizing NOAA-11 AVHRR data, a 
two dimensional scatter plot analysis technique was employed 
to determine and classify the radiative signatures of the 
smoke and dust. A two dimensional mask routine was then used 
to assess the reliability of the scatter plot analyses and 
spatially display the results. A channel 1 to channel 2 
radiance ratio and a channel 4 brightness temperature 
combination provided the best separation of the smoke 
signature from water. The dust plume was unambiguously 
represented by a channel 5 brightness temperature minus 
channel 4 brightness temperature image and a channel 4 
brightness temperature combination. Together the 2D scatter 
plot technique and 2D mask form the groundwork for a possible 
detection algorithm. 
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INTRODUCTION 



I . 



A. MOTIVATION 

On 1 March 1991, the region of the Arabian Peninsula and 
Arabian Gulf presents a unique opportunity to remotely study 
the characteristics of two categories of complex aerosols: 
oil smoke and desert dust. A smoke plume is being produced by 
some 610 burning Kuwaiti oil wells ignited by the retreating 
Iraqi army during the recent Gulf War. Estimates by Johnson 
et al . (1991) put the total burn rate of the oil fires at 

202.5 metric tons per year or 3.9 million barrels per day. 
Additionally, a dust storm on the southern Arabian Peninsula 

created a large dust plume with dimensions several hundreds of 
kilometers in length and width. The presence of these two 
types of aerosols in the same region on the same day offers an 
unparalleled "data rich" environment ideal for a study in 
satellite detection of complex aerosols. 

Both military and scientific concerns motivate the need 
for accurate aerosol detection. First, aerosol concentrations 
can negatively impact military operations. They reduce 
visibilities, affect weapon, machine and human effectiveness 
and hamper aircraft/satellite intelligence gathering efforts 
vital to most military operations. Accurate real-time 
analyses of aerosol spatial distributions and concentrations 
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would be a valuable mission planning aid for all military 
commanders. Second, aerosols can effect changes in local, 
regional and in some cases the global environment. Some of 
these effects include sharp increases and decreases in surface 
and atmospheric heating rates, production of anomalous 
mesoscale weather features and long-term changes in 
climatology. Thus, knowing distributions and concentrations 
of aerosols can aid scientists in studying and predicting 
aerosol effects on the environment. 

B . GOALS 

The primary goal of this project is to accurately identify 
the radiative signatures for smoke and dust over the Arabian 
Gulf and Arabian Peninsula on 1 March 1991. To achieve this 
goal, a scatter plot analysis technique is employed. This 
technique is currently used for cloud classification wherein 
cloud type is identified by a combination of its IR brightness 
temperature and its solar reflectance. The actual detection 
and classification of smoke and dust, however, is complicated 
by the spectral radiative influence of atmospheric 
constituents (primarily water vapor), clouds and geographic 
features (land and water) on the aerosol radiative signatures. 
In particular, it is difficult to detect smoke over water and 
dust over land. Figure 1, a 1 March 1991 visible satellite 
image of the Arabian Peninsula/Gulf, demonstrates these 
difficult effects. Smoke being generated by the oil fires is 
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Figure 
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) . Dust plume is 
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discernible over land but not over the gulf (area 1) and the 
dust storm over the southern Arabian Peninsula (area 2) is 
barely discernible, if at all, from the surrounding desert. 
Therefore, a major challenge in this project is to determine 
the spectral radiance of aerosol particles and then isolate 
them using the scatter plot technique. 

In addition to the primary goal, there are three secondary 
goals embedded in this project. The first is to evaluate the 
scatter plot analysis technique as to its viability for 
analyzing aerosol radiative signatures. Second, the project 
will evaluate the two dimensional masking routine as to its 
viability for representing and analyzing aerosol spatial 
distributions and concentrations. And third, these aerosol 
signatures will aid in development of an aerosol detection 
algorithm for possible military application. 
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II. BACKGROUND 



A. BASIC THEORY 

The theory behind this study is based on the fundamentals 
of radiative transfer. Simply by knowing what can happen to 
a photon in a medium (as a function of its wavelength) and 
what photons a satellite detects, one can understand the basis 
of the analyses. 

A photon will either be absorbed by the medium, reflected 
by the medium or transmitted through the medium. Again, what 
happens depends on the wavelength of the photon and the 
physical characteristics of the medium. A satellite will 
detect photons reflected by the surface and/or aerosols or 
those emitted by the surface, aerosol and/or atmosphere. In 
either case, the number of photons detected is dependent upon 
the physical characteristics of the surface, aerosol and 
atmosphere as a function of wavelength (JL). 

The following simple equation summarizes the satellite 
detection where L = monochromatic radiance: 

~ Ii«uirCAC*(X) Lia«ito«ox(X) 

Aerosols are detectable when the aerosol radiance contrast 
strongly with the ambient surface and/or atmospheric radiance. 
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B. CHARACTERISTICS OF OIL SMOKE 



The Kuwaiti oil fires represent the first time a large 
scale oil smoke plume has been available for scientific 
research. Therefore, little is known about the actual 
radiative effects and characteristics of a smoke plume of this 
type. Laboratory experiments, however, have shown 

approximately 100% of aerosol from burned oil is elemental 
carbon (Crutzen et al . , 1984). Elemental carbon is highly 
absorbent throughout the solar spectrum which is why oil smoke 
generally appears black at visible wavelengths (Turco et al., 

1990) . Additionally Turco et al. (1990) has noted that smoke 

radiative absorption decreases slowly with increasing 

wavelength . 

Actual measurements of the Kuwaiti smoke plume, obtained 
by a British Meteorological Office Research flight on 18 March 
1991, indicates near field particles are composed by spherical 
particles of approximately l)im diameter formed into aggregates 
up to several micrometers in length/size (Johnson ,et al., 

1991) . Additionally, there is high near source water vapor 

content giving the potential for condensed water drops. The 

water vapor is a combustion byproduct as a result of water 

intrusion into the oil wells (Limaye et al . , 1991). Finally, 
there are reports of oil drops in the near field plume. 
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C. CHARACTERISTICS OF DESERT DUST 

A study by Carlson and Benjamin (1980) identified some 
general characteristics of dust. First, it is highly 
reflective through all visible wavelengths. Second, there is 
a general trend of increasing absorption with increasing 
wavelength from near IR into middle IR. This is indicated by 
Figure 2 which gives a distribution of the complex index of 
refraction for desert aerosols. Durkee (1984) provides a good 
discussion of how the imaginary index of refraction controls 
the absorption by particles. Finally there is a high 
variability in absorption in the IR water vapor window (8-12 
microns). Figure 3, which is an enlargement of Figure 2, 
demonstrates this characteristic. Notice the variability 
between the channel 4 and 5 wavelength bands. Dust absorbs 
channel 4 wavelengths better than channel 5. This is an 
important characteristic for dust detection in this study. 

D. 1 MARCH 1991 SMOKE/DUST PLUME DISTRIBUTION 
1. Smoke Plume 

There are two layers of smoke as a result of the 
Kuwaiti oil fires. A low level smoke plume is being advected 
southeastward over the Gulf and along the Saudi Arabian coast. 
In the vicinity of Qatar, the plume is turned and advected 
inland over Saudi Arabia. This was determined by analyses of 
various AVHRR passes along with surface and upper air wind 
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and pressure fields. The weather data was obtained from Fleet 
Numerical Oceanography Center (FNOC). Figure 4 shows the 0900 
UTC surface pressure analysis and surface wind observations 
which support the advection pattern. The top of the low level 
plume is between 1650 and 2115 meters. This was crudely 
determined utilizing the 1 March 1991 00 UTC Dharhan, Saudi 
Arabia sounding. Figure 5, (location on Dharhan provided later 
in Figure 13) along with a 1 March 1991, 0352 UTC channel 4 
brightness temperature image. Essentially, a temperature 
range for the plume tops was extracted from the brightness 
temperature images and plotted on the sounding temperature 
profile. The plotted temperature correspond to a particular 
height in millibars which was converted to meters. The 
estimated altitude range falls within the observed smoke layer 
range of 1000 to 3000 meters as reported by Limaye et al . 
(1991). The data for the Limaye et al. study was taken on 18 
March 1991 under similar meteorological conditions as on 1 
March . 

An upper level smoke layer between 3600 and 4200 meters is 
advected eastward over the Gulf and towards Iran. This 
correlates well with strong westerly winds 700 millibars and 
above as indicated in Figure 5. The altitude range of this 
upper level smoke plume was not determined in the project, but 
comes from the reported observations in the Limaye et al. 
(1991) study. 
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Figure 4 . 1 March 1991 0900 UTC surface pressure analysis in 

millibars and station surface winds in knots. 



PRESSURE (mb) 




Figure 5. 1 March 1991 00 UTC Dharhan, Saudi Arabia sounding. 
Used to determine low level Ismoke plume between 830mb to 
785mb or approximately 1650m to 2115m. 
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2. Dust Plume 



The dust plume has its source approximately 335km 
south southeast of Riyadh, Saudi Arabia. The plume is 
adverted southwestward approximately 515km to the base of the 
Hejaz Asir mountains. This advection is supported by the 
surface pressure pattern seen in Figure 4 . The dust cloud 
then spreads some 600-700km along the base of the mountains. 
Tops of the dust plume range from 3000 to 4000 meters. The 
tops were determined in the same way as the low level smoke 
plume except the 1 March 1991 12 UTC Riyadh, Saudi Arabia 
sounding. Figure 6, and a 1 March 1991 1026 UTC channel 4 
brightness temperature image were used. 
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Figure 6. 1 March 1991 12 UTC Riyadh, Saudi Arabia sounding. 
Used to determine dust plume top between 700mb and 650mb or 
approximately 3000m to 4000m. 
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III. PROCEDURES 



A. DATA 

Satellite data used in this project was collected by the 
National Oceanographic and Atmospheric Administration (NOAA) 
10 and 11 sun synchronous polar orbiting satellites. The NOAA 
10 data was taken on 1 March 1991 at 0352 UTC (0652 local) in 
a descending node over the Middle East and the NOAA 11 data 
was taken on 1 March 1991 at 1026 UTC (1426 local) in an 
ascending node over the Middle East. Both the NOAA 10 and 11 
data are contained on National Environmental Satellite Data 
and Information Service (NESDIS) tapes obtained from the Navy 
Oceanographic and Atmospheric Research Lab West (NOARL West). 
Advanced Very High Resolution Radiometer (AVHRRs) measured the 
upwelled radiance. AVHRR/1 on NOAA 10 has four channels or 
spectral bands and the AVHRR/2 on the NOAA 11 satellite has 5 
channels. Table 1 indicates the bandwidth of the two 
radiometers . 

NOAA 10 data was used to determine the approximate 
vertical extent of the low level smoke plume. NOAA 11 data 
was used to derive parameters for all scatter plots and to 
determine the vertical extent of the dust plume. 
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Table 1. AVHRR/1 AND 2 CHANNEL BANDWIDTHS (ADAPTED FROM 
KIDDER AND VONDER HAAR, 1992) 



CHANNEL 


AVHRR/1 BW 


AVHRR/2 BW 


RADIANCE 


1 


. 55- . 68 )lm 


. 58 - . 684 


Red Visible 


2 


.75-1. 104 m 


. 725-1 . 104 m 


Near IR 


3 


3.55-3. 934 m 


3.55-3. 934 m 


Middle IR 


4 


10 . 5-11. 54 m 


10 . 3-11. 34 m 


Thermal IR 


5 


Channel 4 
Repeated 


11 . 5-12. 54 m 


Thermal IR 



B. SATELLITE DATA PROCESSING 

All processing was performed in the Interactive Digital 
Environmental Analysis Laboratory (IDEA Lab) at the Naval 
Postgraduate School (NPS) in Monterey, California. Initial 
satellite data processing utilized two IDEA Lab programs. 
AVIAN 3.0 written by Motell et al., (1991) was used to glean 
tapes and produce desired satellite parameters/images from 
satellite overviews. And DIPS, a program designed to work 
with existing satellite images, was used to perform 
enhancements and annotate on images. DIPS is essentially a 
"finishing" program for satellite images. The TEKTRONIX RGB 
color printer was utilized for all image reproduction via the 
IDEA Lab workstations. 
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C. DATA ANALYSIS 

1. Radiative Signature Analysis 

Seven parameters/ images are derived from the NOAA 11 
AVHRR/1 to assist in analyzing feature signatures. A two 
dimensional scatter plot analysis technique is then employed 
to analyze the individual signatures. After the feature 
signatures are determined, a two dimensional mask program is 
used to spatially display the results of the analyses and 
assess their reliability. 

a. Derived Parameters 

The following seven parameters are utilized in this 

project : 

1. Channel 2 calibrated albedo (Cal 2) 

2. Channel 1/Channel 2 radiance ratio (S12L) 

3. Channel 3 Brightness Temperature (Temp 3) 

4. Channel 4 Brightness Temperature (Temp 4) 

5. Channel 5 Brightness Temperature (Temp 5) 

6. Temp 3 minus Temp 4 (T3-T4) 

7. Temp 4 minus Temp 5 (T4-T5) 

From this point forward, the parameters will be referred to by 
their abbreviations in parenthesis. 

Cal 2 provides an indication of the solar 
reflectance of image features. S12L is used to determine 
smoke distributions over water. It uses the differing 
absorption character ist ic of the Gulf water between channels 
1 and 2 to contrast with the smoke signature. Temp 3, during 
the daytime, contains both reflected solar energy and emitted 
IR. Temp 4 wavelength band is in a water vapor window and 
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therefore gives a relatively accurate indication of feature 
temperatures. Temp 5 also gives feature brightness 
temperature. Its wavelength band, however, is located just on 
the edge of the water vapor window so there is some 
sensitivity to water vapor in the feature brightness 
temperatures. T3-T4 gives a first order estimate of the 
reflected solar radiation in Temp 3. T4-T5 has been typically 
used to sense water vapor effects on sea surface temperature 
measurements. In this project, however, it is used to 
separate the dust from other features based on the differing 
dust absorption characteristic between channels 4 and 5 as 
previously discussed in Chapter II. 

b. Scatter Plot Technique 

In the production of a full image, two dimensional 
scatter plot is a multi-step process. Figure 7 outlines the 
steps in the procedure. First, two satellite subimages or 
parameters are selected. Subimages are 512 by 512 picture 
element (pixel) areas derived from larger overview regions of 
a satellite pass. The data for each pixel in the subimages is 
stored in ascii format and must be converted to real numbers. 
This is done by a Fortran program called "Image to Real". The 
program converts the ascii data for an image into a 512 by 512 
real number array data file. Once the conversion is complete, 
both arrays are input into a Fortran program called "Get 
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Figure 7. Procedural flowchart for producing full image 
2D scatter plot for two subimages. 
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Data", which plots image one versus image two on a 2D scatter 
plot. Note, the number of pixels plotted from each image can 
be varied. Figure 8 is an example of a full image 2D scatter 
plot with every other pixel plotted. 

The next procedure, after producing a full image 2D 
scatter plot, is to identify the individual feature radiative 
signatures on the plot. The flow chart (Figure 9) outlines 
this process. Utilizing a Fortran program called "Store", 
individual feature areas (e.g., areas of only land, smoke, 
dust, etc.) can be sampled from a subimage. "Store" then puts 
the ascii information from the area into a data file. That 
process is repeated until enough areas are sampled to 
accurately identify the features of interest on the full 2D 
scatter plot. Once the desired sample area ascii data files 
are created for each subimage (keeping the sample area 
locations consistent for each) individual sample area scatter 
plots are produced. Figure 10 is an example of a sample area 
scatter plot. This is done by inputting the individual sample 
data files into the Fortran program "Get Ascii". "Get Ascii" 
produces the same scatter plots as "Get Data" except it plots 
only the sample areas from each subimage versus each other. 
By doing this, the general location and approximate 2D 
radiative signature of a desired feature can be analyzed on 
the full image scatter plot. Each sample area is 
approximately represented by either an oval or a hand drawn 
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Figure 8. Example of full image 2D scatter plot of 
versus Temp 4 for area 1 in Figure 14. 
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Figure 9. Procedural flow chart for producing image sample 
areas on 2D scatter plot. 
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CH4 TEMP. (K) SMOKEl (iO.S-fxU) 



Figure 10. Example of single feature 2D scatter plot. Smoke 
1 corresponds to SI sample area on Upper Gulf Images, see 
Figure 16. 
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sketch on the full image 2D scatter plot. The goal is to 
encompass at least 80 percent of the pixels defining a sample 
area. A hand drawn sketch is used when either the shape of 
the pixel area representing a feature is not conducive to an 
oval and/or if 80 percent of the feature cannot be encompassed 
without highly exaggerating the pixel area defining the 
feature. The drawing of the ovals and/or hand sketches is 
subjective and in most cases will not exactly align with the 
actual analyzed signature limits of the features. They are 
drawn to give an idea how the feature signatures breakout on 
the full image scatter plot. 

c. Verification of Radiative Signature Analyses 

To verify the signatures analyzed using the scatter 
plot technique, a two dimensional mask program called "Maskit" 
is utilized. Figure 11 is a procedural flow chart showing 
this process. Essentially, the analyzed two dimensional 
limits of a feature are input into "Maskit". The program then 
colors or masks any pixels falling in the limits with one of 
seven primary colors. The masked area provides the feedback 
to validate the analyzed results. Figure 12 is an example of 
analyzed limits on a scatter plot and Figure 13 is the mask of 
those limits. 
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Figure 11. Procedural flow chart for verification 
of radiative signature analyses. 
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IMAGE CLASSIFICATION 




Figure 12. Same as Figure 8 except example of analyzed 2D 
signature limits on full image 2D scatter plot of S12L and 
Temp 4 for area 1 in Figure 14. Region represented in this 
image corresponds with area 1 of Figure 13. 
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Figure 13. Example of mask 
12. Red is thick smoke. 
Magenta is thin smoke over 



using analyzed limits from Figure 
Green is thin smoke over land, 
water. Blue is water. 
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IV. SCATTER PLOT ANALYSES/RESULTS 



A. UPPER GULP SCATTER PLOTS 

The images, used to produce the scatterplots for analysis, 
cover the same area 1 of the upper Arabian Gulf (see boxed 
area on map Figure 14). Each image is 512 by 512 picture 
elements (pixels) in size with a l.llcm per pixel resolution. 
The images were generated from a NOAA 11 AVHRR Middle East 
pass ta)cen at 1026 UTC 1 March 1991. The actual scatter 
plots, although taken from the 512 by 512 images, were 
generated by sampling every other pixel. Thus the scatter 
plots effectively have a 2km resolution. Sampling every pixel 
proved to be too cumbersome and made it generally more 
difficult to discern the desired signatures. 

The following features were analyzed on each scatter plot 
to determine their radiative signature: 1) water, 2) thick 
smoke, 3) thin smoke over land, 4) thin smoke over water and 
5) land. The "water" signature represents the portions of the 
gulf without a discernible smoke or cloud overcast. Those 
areas of "clear" water are primarily concentrated in the far 
northern gulf and along the Iranian coast. "Thick smoke" 
refers to smoke with a radiative signature the same as smoke 
within approximately 10 to 100km of the main oil well fires. 
As might be expected, the thick smoke areas stretch from the 
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sources down the central portions of the various plumes. The 
"thin smoke over land" signature is representative of the land 
areas lust perceptible through the smoke overcast. Those 
areas are concentrated on the fringes of the smoke plumes. 
The signature for "thin smoke over water" depicts any water 
area with a "thin" smoke overcast. The thin smoke over water 
primarily extends from the plume fringes in the northwest 
Arabian Gulf through the central gulf down to the southeast 
where the signature is lost due to dispersion. "Land" with no 
smoke or cloud overcast accounts for most of the radiative 
signature indicated on the scatter plots. As will be noted in 
the following analyses there are essentially four discernible 
land signatures. Figure 15 presents AVHRR channel 1 to 
channel 2 radiance ratio image (S12L) of the upper gulf and 
depicts, with colored boxes, the sampled pixel areas of the 
features previously listed and described. The S12L and Temp 
4 (Figure 16) images served as the "baseline" images to 
determine feature sample areas. Note some features were 
sampled in more than one location. As a reminder, on all 
scatter plots, ovals and/or hand drawn sketches encompass the 
major itv of those pixels which are representative of a 

particular feature. They are not meant to represent the exact 
analyzed signature limits. Finally, there was experimentation 

with a number of image combinations. Only the most 
illuminating cases, however, are presented. 



30 




Figure 15. 1 March 1991 1026 UTC Upper Gulf S12L image 
depicting feature sample areas (L-land, S-smoke, W-water, TSW- 
thin smoke over water and TSL-thin smoke over land). 
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as in Figure 15. 
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1- Case 1 Scatter Plot Analysis; S12L Versus Temp 4 
(Figure 17) 
a. Water: 

Water ranges from approximately 3.06 to 5.0 on the 
S12L axis and from 287.29°K to 292. 5°K on the channel 4 
temperature axis. "Water" on the S12L image constitutes the 
brightest feature (mid gray to white). This implies the gulf 
water reflects more radiation in channel 1 than in channel 2. 
At this time it should be noted, theoretically given "pure" 
water with a relatively deep depth (>~10m) the S12L ratio 
would be close to one. This results from the high absorption 
for both channels 1 and 2 wavelength bands and because 
reflectance is low and approximately the same. There are two 
major factors which account for the differences in reflectance 
and thus give a relatively high range of radiance ratios. 
First is the amount of sediment in the water. February and 
March are climatologically in the middle of the rainy season 
for the Middle East as indicated (Taha et al . , 1981). 
Consequently, it is assumed there is sediment laden river 
runoff, primarily in the northern gulf around the Euphrates 
and Tigris River deltas and in the northwestern gulf, along 
the coast of Iran, due to various rivers flowing down from the 
Zagros Mountains. The large amount of sediment in those 
coastal waters increases the reflectance of the water in both 
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Figure 17. Case 1 Upper Gulf scatter plot. Enclosed areas in 
Figures 15 and 16 (W-water, L-land, S-smoke, TSW-thin smoke 
over water and TSL-thin smoke over land). 
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channels 1 and 2. But because channel 2 has a wavelength band 
at the end of the visible spectrum and into near IR, the water 
absorbs more at those wavelengths compared to channel 1, which 
has a wavelength band wholly in the visible spectrum. In 
fact, channel 1 wavelengths are absorbed approximately 5 times 
less than those of channel 2. 

The second factor accounting for the radiance 
variation between channels 1 and 2 is the bathymetry of the 
gulf. Although the variation is not as dramatic as that 
attributed to the sediment, there is a correlation between 
water depth and the magnitude of the S12L ratio. As the water 
depth increases the radiance ratio decreases. Especially 
evident of this correlation is the outline of the 10 fathom 
line on the S12L image. The explanation for this correlation 
appears to be twofold. First, near surface suspended sediment 
correlation decreases as a function of distance from shore due 
to particle settling and dispersion. Since the amount of 
suspended sediment seemingly drives the amount of channel 2 
waveband absorption (i.e., the more sediment the more 
absorption) it ma)<es sense that channel 2 and channel 1 
radiances would get closer, meaning the magnitude of the S12L 
ratio would decrease - as it does. Second, with increasing 
water depth, there is increased absorption in both channels. 
In other words, once a photon from either waveband enters the 
water, the deeper the water the less chance the photon has to 
be scattered bac)< out. Therefore, as water depth Increases, 
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the observed radiances in channels 1 and 2 decrease and begin 
to equalize. This translates to lower and lower S12L ratios 
or darker and darker gray shades on the S12L image. The 
channel 4 temp range for water, as indicated earlier, is 
approximately 286°K to 292. 5°K. As might be expected the 
cooler water is in the northern gulf and warmer water in the 
South. Latitudinal air temperature variations as well as the 
influx of cool river water in the northern gulf are suspected 
as the primary reasons for the 6.5°K temperature difference. 
i>. Thick Smoke: 

Thick smoke ranges from approximately 2.10 to 2.37 
on the S12L axis and from 282. 2°K to 288.11°K on the channel 
4 temperature axis. The S12L ratio of around ”2" implies the 
smoke particles are approximately twice as reflective of the 
channel 1 wavelength band compared to channel 2's. That, 
however, is most likely not the case. Due to intrusion of 
water into the oil wells, water vapor, in addition to smoke 
particles and a myriad of other gases, became a combustion 
product (Limaye et al., 1991). As a result, the smoke 
particles act as condensation nuclei and introduce liquid 
water into the smoke plume. Since water is slightly more 
absorptive of channel 2 wavelengths, observed channel 1 
radiance is slightly higher which accounts for an S12L ratio 
greater than 1 . 
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The channel 4 thick smoke temperatures should be 
very representative of the actual temperatures at the top of 
the smoke layer. USS Wisconsin (BB 64) (1991), post cruise 
report indicates total obscuring of the sun near the source(s) 
which likely translates to very little IR contamination/ 
contribution to the observed smoke top temperature due to 
surface/ground emittance. 

c. Thin Smoke Over Land: 

Thin smoke over land ranges from approximately 1.61 
to 2.09 on the S12L axis and from 282. 6®K to 292.59°K on the 
channel 4 temperature axis. The lower S12L ratios (compared 
to thick smoke) are primarily due to the surface/ground 
reflectance contribution to the thin smoke signature. The 
observed surface/ground radiances are approximately the same 
in channels 1 and 2 which implies an S12L ratio of around 1. 
Thus, the "ground contamination" decreases the thin smoke 
ratio . 

Surface/ground IR emittance also accounts for the 
slightly higher channel 4 temperature range for the thin 
smoke. Compared to totally obscured, as near the sources, the 
ground receives some solar radiation and is heated slightly. 
The heating results in surface IR emittance and a larger 
temperature range for the thin smoke over land. 
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d. Thin Smoke Over Water: 



Thin smoke over water ranges from approximately 
2.38 to 3.05 on the S12L axis and from 285.25°K to 291.17°K 
on the channel 4 temperature axis. The signature for thin 
smoke over water is strongly influenced by the water over 
which it lies. As in the case of thin smoke over land the 
water radiance contaminates the basic smoke signature and this 
accounts for the higher S12L ratio and warmer temperature 4 
range - as compared to thick smoke. 
e. Land: 

Land primarily ranges from 1.0 to slightly above 
2.5 on the S12L axis and is spread across the full range of 
channel 4 temperatures with most of the signature concentrated 
above 290°K. The reason for the diversity in signature rests 
in the fact there are essentially four different land features 
in the images. The features vary from desert, 

irrigated/vegetated plots, river deltas and finally snow 
covered mountainous areas. In general, however, the land 
areas are the warmest features in channel 4 temp and the 
darkest (lowest ratio) in S12L. The warm signature occurs 
because the satellite data was taken at 1426 local time - the 
land should be close to its warmest diurnal temperature. The 
relatively small range in S12L implies most land features have 
the same reflectance in channel 1 as in channel 2. 
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f. Sutamary: 

The S12L versus channel 4 temp scatter plot 
highlights the signatures of the five main features quite 
well. Water is confined to a relatively narrow low to mid 
axis temperature range of 6.5°K but accounts for the highest 
S12L ratios. This translates to water appearing as a light 
gray on the temperature 4 image and mid level gray to white on 
the S12L image. Thick smoke lies at the low end of the 
temperature axis and covers about a 6°K temperature span. Its 
S12L image is very narrow (~.3) and is centered around a ratio 
of 2.23. Thick smoke appears as a very light to almost white 
gray shade on the temperature 4 image and as dark gray on the 
S12L image. Thin smoke over land has a temperature signal 
spread across the lower half of the temperature 4 axis but has 
a relatively narrow (~.5) S12L range centered at 1.85. Its 
overall signature is driven by the land it covers, which 
accounts for the wide temperature range and narrow S12L ratio. 
It appears across the full spectrum of mid-level grays in the 
temperature 4 image and as a dark gray (not black) on the S12L 
image. Thin smoke over water covers a 6°K, low to mid, 
temperature range and lies in a . 7 S12L ratio range centered 
at approximately 2.73. As with its counterpart over land, the 
water has a significant influence on its signature. Thin 
smoke over water appears as a light gray on the temperature 4 
image (almost indiscernible from the water) and as a mid level 
dark gray on the S12L image. The highly variable land 
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signature is also the most dominant feature on the scatter 
plot. It lies across the entire temperature axis and accounts 
for the lowest S12L ratios. It appears anywhere from the 
darkest gray (black is the actual oil fires) to white on the 
temp 4 image and from black to dark gray on the S12L image. 

Overall, this scatter plot separates the signatures 
of the five main features the best. Its most important 
aspect, however, is the separation of thin smoke over water. 
Even though the range limits are highly subjective, this is 
the only plot that discernibly separates the thin smoke over 
water from the other features. 

2. Case 2 Scatter Plot Analysis: Cal 2 Versus Temp 4 

(Figure 18) 
a. Water: 

Water ranges from approximately .7 to 1.7 percent 
albedo on the channel 2 axis and from 287.29°K to 292. 5°K on 
the channel 4 axis. The channel 4 temperature range was 
previously discussed in case 1. The channel 2 albedo range 
confirms what is generally known about water - that it has a 
low albedo especially when the sun is close to its zenith, as 
in this case (Kidder and Vonder Haar, 1992). Therefore, water 
appears black in the channel 2 image. 
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Figure 18. Case 2 Upper Gulf scatter plot. Feature 
designation as in Figure 17. 
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t>. Thick Smoke: 



Thick smoke ranges from approximately 1.7 to 2.9 
percent albedo on the channel 2 axis and from 282. 3°K 
to 288.11°K on the channel 4 axis. The channel 4 range is 
consistent with that analyzed in case 1. The low albedo of 
smoke in channel 2 implies the smoke particles are highly 
absorbent of incident solar radiation. In fact, almost no 
visible light is reflected from the particles. The light that 
is reflected is most likely a result of scattering by the 
condensed water on the larger smoke particles. 

c. Thin Smoke Over Land: 

Thin smoke over land ranges from approximately 2.0 
to 4.0 percent albedo on the channel 2 axis and from 282. 9°K 
to 287. 8°K on the channel 4 axis. The albedo signature is 
higher than the thick smoke albedo because the particle 
density is less which allows visible light reflected from the 
primarily sand surface to boost the albedo slightly. The 
temperature range is consistent with that analyzed in case 1. 

d. Thin Smoke Over Water: 

Thin smoke over water signature ranges from 
approximately 1.4 to 2.1 percent albedo on the channel 2 axis 
and from 285.25“K to 291°K on the channel 4 axis. Note, these 
ranges are the most subjective of all the features. Because 
of the low albedo of both smoke and water it is extremely 
difficult to discern a cutoff. Considering the resolution of 
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the images and scatter plots, the cutoff is probably 
indiscernible. Overall, however, the albedo signature is 
generally lower than both the thick smoke and thin smoke over 
land signatures. This is due to the influence of the very low 
water albedo. The temperature 4 signature is consistent with 
that analyzed in case 1. 

e. Land: 

The land signature ranges from approximately 2.0 to 
30.0 percent albedo on the channel 2 axis and covers the full 
temperature range on the channel 4 axis. There were some 
pixels higher than 30 percent albedo but the channel 2 axis 
was limited to 30 to enhance the "main" image features. Those 
pixels higher than 30 percent albedo constitute areas of 
highly reflective snow, in the Zagros mountains (upper right 
corner of image. Figure 19), and some small pastches of desert 
sand on the west coast of the gulf (lower left corner of image 
Figure 19). The lower albedos represent the 
vegetated/irrigated areas. Once again, the temperature 4 
signature is consistent with that analyzed in case 1. 

f. Summary: 

Overall, the scatter plot indicates smoke and water 
have similar albedo characteristics. Water, thick smoke and 
thin smoke over water/land are concentrated between .7 to 4.0 
percent albedo. This translates to very little separation. 
Water is at the low end of the albedo range and thin smoke 



43 




Figure 19. 
Br ightness 
areas with 



1 March 1991 1026 UTC Upper Gulf Cal 2 image, 
pixels in upper right and lower left corners are 
an albedo greater than 30 percent. 
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over land is at the high end. The temperature differences of 
the features account for most of the small signature 
separation observed. The "standout" signature on this plot is 
that of the land. In fact, there is excellent separation of 
the individual land regions. This separation especially 
emphasizes the diversity of land in the upper gulf region. 

3. Case 3 Scatter Plot Analysis: T3-T4 Versus Temp 4 

(Figure 20) 

In channel 3 there is contribution to total daytime 
radiance from both the solar and thermal IR spectrums. Thus, 
channel 3 daytime brightness temperatures do not represent 
just thermal emission but are "contaminated" with some 
reflected solar radiance. For this reason, observed channel 
3 daytime brightness temperatures are generally higher than 
those observed in a strictly thermal IR channel such as 
channel 4. Therefore, when a channel 3-4 temperature 
difference image is produced most values are positive as can 
be seen on the scatter plot, 
a. Water: 

Water ranges from approximately -.8°K to 4.0°K on 
the T3-T4 axis and from 287.29°K to 292. 5°K on the channel 4 

temperature axis. It has the lowest difference range of the 
main features. This is consistent with the low albedo of 

water. Because water does not reflect much solar radiation 
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Figure 20. Case 3 Upper Gulf scatter plot. Feature 
designation as in Figure 17. 
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its brightness temperature in channel 3 are close to those in 
channel 4. The upper limit is most likely indicative of the 
shallower, sediment laden water and the lower limit indicative 
of the deeper, reduced sediment water due to the differences 
in reflectance as discussed in case 1. The channel 4 
temperature range is consistent with those analyzed cases 1 
and 2. The temperature bounds, however, were more definitive 
in this scatter plot and thus easier to analyze, 
f). Thick Smoke: 

Thick smoke ranges from approximately 2.2°K to 
9.0°K on the T3-T4 axis and from 282. 2°K to 288.11°K on the 
channel 4 axis. The bottom limit on T3-T4 axis is relatively 
discrete but the upper limit is highly subjective which 
emphasizes the complexity of the thick smoke areas. Due to 
the reflective component in channel 3, the thick smoke 
signature is dispersed over a wide range of T3-T4 values. 

The thick smoke sample areas (see SI and S2 boxes 
on Figure 21) encompasses smoke over both land and water and 
part of the "fringe" areas of the smoke plumes. The land 
(sand) in that area has a relatively high albedo. That 
reflective influence in the fringe areas furthest from the 
fires, is what drives the channel 3 temperature higher and is 
responsible for the higher T3-T4 values at the warmer channel 
4 temperatures (286°K - 287°K range). The higher T3-T4 values 
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21. 1 March 1991 1026 UTC T3-T4 
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Upper Gulf image, 
land and water . 
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in the cooler channel 4 temperature range (283°K - 284°K) are 



probably driven in a small part by the land reflectivity but 
in a large part due to the back scattered visible radiation by 
the water drops attached to the smoke (this is in the area 
nearest the fires where the atmospheric water vapor content is 
the highest). Note this is purely speculation to account for 
the two "thick smoke branches" seen on the scatter plot. The 
thick smoke values in the 3.8°K to 5.6°K range on the T3-T4 
axis can most likely be accounted for by the varying amounts 
(albeit small) of added reflectance in channel 3 from the gulf 
water . 

The lowest T3-T4 values represent the thickest 
smoke. The denser and thicker the smoke, the less chance 
there is for reflective influence in channel 3 from the land 
or water which implies the 3-4 temperature difference should 
be low. The thickest smoke channel 4 temperature spread is 
primarily due to whether the smoke lies over land (cooler 
temps) or water (warmer temps). The channel 4 temperature 
range is consistent with that analyzed in the previous cases, 
c. Thin Smoke Over Land: 

Thin smoke over land ranges from approximately 
2.0°K to 16.0°K on the T3-T4 axis and from 282. 9°K to 289. 8°K 
on the channel 4 axis. As with thick smoke, it is dispersed 
over a wide T3-T4 range. The reason is twofold. First, as 
discussed under thick smoke, the high albedo of the sand 
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surface contributes to the anomalously high channel 3 
brightness temperatures which cause the positive T3-T4 
difference values. Second, the T3-T4 range is large because 
the contribution by surface reflectiveness is inversely 
proportional to smoke thickness. The thicker the smoke, the 
more sunlight is blocked from the surface and the lower the 
surface reflectance. There is a possibility that some of the 
lowest T3-T4 values may result not only because of the 
previous relationship but because there may be a large patch 
of vegetated land (farmland) under the thin smoke area. These 
patches would have very low albedos and thus little 
reflectance contribution to channel 3 brightness temperature. 
Note, the actual existence of those vegetated areas is 
unverified. The closeness of these patches to populated areas 
is the cause for speculation. The channel 4 temperature range 
is consistent with the previously analyzed cases. 
d. Thin Smoke Over Water: 

Thin smoke over water ranges from approximately 
3.2°K to 6.0°K on the T3-T4 axis and from 285.25°K to 291.17'»K 
on the channel 4 axis. The influence of the slight reflective 
component of the water is responsible for the larger channel 

3 temperature and thus the positive T3-T4 values. The channel 

4 temperature limits are consistent with the previous case and 
as with case 1 and 2 the limits are very subjective. 
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e . Land : 



The land signature ranges from approximately 1.8°K 
to 26®K on the T3-T4 axis and across the full temperature 
range on the channel 4 axis. Once again, as in case 2, the 
diversity of the land features is illustrated quite well. The 
reflectance of a particular land feature is the key to 
determining its T3-T4 value. As might be expected, the sand 
areas and snow covered mountain areas dominate upper part of 
the T3-T4 axis. Whereas, the river deltas and vegetated areas 
(like L5) are found at the low end of the T3-T4 range. The 
channel 4 temperature range is consistent with cases 1 and 2. 
f. Summary 

The signature separation is indiscrete between 
thick smoke and thin smoke over land/water. There is, 
however, fairly definitive separation of those features from 
the water and land signatures. Water dominates the low T3-T4 
values and has a discrete lower channel 4 temperature limit. 
This is in contrast to the previous two scatter plots (Figures 
17 and 18). Understand, however, this lower limit represents 
gulf water without a smoke overcast. It is expected the water 
overlain by the smoke has slightly cooler temperatures. The 
smoke signature, in general, lies to the far left (lower 
temperatures) on the channel 4 temperature axis but is spread 
across a large T3-T4 difference range of 14°K. It is 
difficult to discern a cutoff between thick and thin smoke. 
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But essentially, the thick smoke is at the lower T3-T4 values 
and the thin smoke at the higher values. The thin smoke 
signatures, over land and water, appear to be strongly 
influenced by surface reflectance. The predominance of the 
sand surface drives the thin smoke overland signature to the 
high T3-T4 values whereas the low albedo of the water brings 
the thin smoke over water signature to the low to mid range 
values ( 3-6 ) . 

Overall, this scatter plot emphasizes the 
significance of the reflective component on the brightness 
temperature in channel 3. The reflectance factor tends to 
make the brightness temperatures anomalously high; as 
indicated by most features on the scatter plot having a 
positive T3-T4 signature. Generally, any feature with a low 
albedo, such as water or smoke falls toward the low end of the 
axis and any with a high albedo, such as snow or sand, falls 
at the upper end. The "standout" signatures as in case 2, are 
the land features. This scatter plot once again accentuates 
the diversity of the land in the upper gulf region. 

4. Case 4 Scatter Plot Analysis: T3-T4 versus Cal 2 

(Figure 22) 

As analyzed in case 3, the influence of the 
reflectance component on channel 3 brightness temperatures is 
what seems to control the T3-T4 difference values. The idea 
of this scatter plot is to compare that reflectance component 
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Figure 22. Case 4 Upper Gulf scatter plot. Feature 
designation as in Figure 17. 
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against the albedos of the features to draw out the actual 
contribution reflectance makes towards increasing the channel 
3 brightness temperatures for each feature. 

a. Water: 

Water ranges from approximately -.8°K to 4.0®K on 
the T3-T4 axis and from .9 to 1.7 percent albedo on the Cal 2 
axis. These ranges are consistent with those analyzed in 
cases 3 and 2. Water has a low albedo and thus its 
reflectance does not significantly increase its channel 3 
brightness temperature. Therefore, the water signature 
resides in the lower left corner of the scatter plot at a low 
T3-T4 difference and a low albedo. 

b. Thi ck Smok e : 

Thick smoke ranges from approximately 2.2®K to 
9.0°K on the T3-T4 axis and from 1.7 to 2.9 percent albedo, 
consistent with cases 3 and 2. The relationship here, 
however, hints the reflectance component is not a major factor 
in increasing the T3-T4 difference for the smoke. There is 
some slight dependence, but comparing the ranges shows a 6.8°K 
T3-T4 difference range compared to a 1 percent albedo range. 
In other words, as the T3-T4 value goes up, albedo changes 
very little - there is not a linear relationship. This 
implies the thick smoke is acting more like a "blackbody" at 
3.7^m (channel 3) than at 10.8jlm (channel 4) because it is 
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emitting more IR radiation at 3.7|im; which corresponds to the 
channel 3 brightness temperatures being higher than in channel 
4 . 

c. Thin Smoke Over Land: 

Thin smoke over land ranges from approximately 
2.0®K to 19.0°K on the T3-T4 axis and from 2.0 to 12.0 percent 
albedo on the channel 2 axis. The upper limits are not 
consistent with case 3 and 2. They are more discernible on 
this scatter plot but are still highly subjective. As with 
thick smoke, the reflectance component does not appear to be 
the key to increasing the T3-T4 difference. The reasons for 
this were discussed under thick smoke. Note, however, the 
thinner the smoke gets, the more linear the relationship 
between albedo and the T3-T4 difference (reflectance 
component). In fact, it becomes almost linear at a T3-T4 
difference of 16®K. At that point the reflectance component 
in channel 3 temperature is the main factor increasing the T3- 
T4 difference. 

Thin Smoke Over Water: 

Thin smoke over water ranges from approximately 
3.2°K to 6.0°K on the T3-T4 axis and from 1.4 to 2.1 on the 
channel 2 axis, consistent with cases 3 and 2. The signature 
is being influenced by the low albedo of the water. But once 
again the characteristic of the smoke particles to radiate 
more IR in channel 3 than channel 4 is controlling the T3-T4 
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difference. The signature tends to blend in with the thick 
smoke signature. 

e. Land: 

The land signature ranges from 1.8°K to 26°K on the 
T3-T4 axis and from 2.0 to 30.0 percent albedo on the channel 
2 axis, consistent with cases 3 and 2. The land demonstrates, 
along with its diversity, the most linearity between the T3-T4 
difference and albedo. The increasing T3-T4 values appear 
highly dependent on the reflectance component in channel 3 
brightness temperatures. Additionally, the signatures of the 
individual land areas are well defined. 

f . Summary: 

The scatter plot has proven to be significant in 
identifying an important characteristic of the smoke 
particles. That is, they seem to absorb/emit channel 3 IR 
wavelengths better than channel 4 wavelengths. This 
characteristic is what controls the increasing T3-T4 values 
for the thick smoke and thin smoke over water signatures and 
the lower T3-T4 range of the thin smoke over land signature. 
The reflective component in the channel 3 brightness 
temperature has very little influence on those signatures. 
Its importance increases, for thin smoke over land, to the 
point where the smoke has dispersed enough to allow the 
surface reflectance to become the key factor in boosting the 
T3-T4 difference. The land and water signatures behave as 
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might be expected. The water is concentrated at the lower 
left corner of the scatter plot (low albedo and low T3-T4 
value) and the land shows a linear correlation between an 
increasing albedo and increasing T3-T4 difference. 



B. ARABIAN PENINSULA SCATTER PLOTS 



The 


Arabian 


Peninsula 


scatter plots were 


produced 


from 


images 


which are 


bound 


by 


the following 


latitudes 


and 


longitudes: 15° 


to 


27° north 


lati^tude by 41 


° to 53° 


east 



longitude (see boxed area on map Figure 14). All images were 
"remapped" from the NOAA 11 AVHRR Middle East overview talten 
at 1026 UTC 1 March 1991. They are 512 by 512 pixels in size. 
But, because they represent a larger area than the Upper Gulf 
images, they have an approximate 2.5km per pixel resolution 
vice a 1.1km per pixel resolution. 

The following seven main features were sampled on each 
image: 1) water, 2) clouds, 3) land, 4) dust, 5) high 
altitude dust, 6) smoke and 7) smoke mixed with and/or 
overlaying dust. Figure 23 shows where these features were 
sampled. Note water and land were sampled in three areas 
each, clouds were sampled in two areas, and smoke, dust, high 
dust and smoke/dust mixed were sampled in only one area each. 
The majority of pixels representing each sample area are 
encompassed on the scatter plots, using ovals and/or hand 
drawn sketches as in the Upper Gulf cases. 
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Figure 23. 1 March 1991 1026 UTC Temp 4 Arabian Peninsula 
image depicting feature sample areas (W-water, L-land, C- 
cloud, D-dust, HD-high dust, S-smoke, SDmi xed-smok e and dust 
mixed ) . 
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The primary goal of this paper, as stated in the 
introduction, is to identify the radiative signatures of smoke 
and dust. Therefore, only dust, smoke, high dust and 
smoke/dust mixed are strictly analyzed in the following 
scatter plot analyses. Much of the Upper Gulf discussion 
concerning the land and water radiative signatures is 
applicable here, however. 

For clarification, "dust” refers to low altitude dust and 
"high dust" refers to high altitude dust. Additionally, the 
"smoke" area analyzed on the Arabian Peninsula images is 
similar in location to the "thin smoke over land" area 
analyzed in the Upper Gulf images. A number of image 
combinations (cases) were generated. The following cases are 
the most illuminating. 

1. Case 1 Scatter Plot Analysis: Cal 2 versus Temp 3 

(Figure 24) 

This image combination was selected because it 
highlights the significant difference in solar reflectance 
between smoke and dust. 

3 - Dus t : 

Dust ranges from 25.5 to 30.2 percent albedo on the 
Cal 2 axis and from 320.8° K to 321.8° K on the Temp 3 axis. 
The relatively high albedo is consistent with the general dust 
characteristics previously discussed in Chapter II. The 
temperatures in channel 3, however, are anomalously high. In 
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Figure 24. Case 1 Arabian Peninsula scatter plot, 
designation corresponds with sample areas in Figure 



Feature 

23. 
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Carlson and Benjamin's (1980) idealized case for Saharan dust, 
it can be inferred that dust at 3.7)im has a high optical 
depth, small imaginary index of refraction (~.01), and high 
single scatter albedo. These characteristics agree with the 
channel 3 temperature signature in this case. The dust 
temperature is being driven abnormally high by the large 
reflective component unaccounted for in the brightness 
temperature. It is highly unlikely the temperature at the top 
of the dust plume is 39°C. 
b. High Dust: 

High dust ranges from approximately 26.5 to 34.0 
percent albedo on the Cal 2 axis and from 320. 6°K to 321. 8°K 
on the Temp 3 axis. The slightly higher albedo, compared to 
that of lower altitude "dust", may be attributed to two 
factors. First, it is reasonable to assume the high dust 
consists of generally smaller particles than lower level dust. 
Consistent with Mie scattering theory, decreasing particle 
size while keeping other parameters relatively constant has 
the effect of increasing upwelled radiation (Kidder and Vonder 
Haar , 1991). Second, the area of the high dust is several 

hundred kilometers down wind of the plume source. Thus, it is 
likely there is more wind generated aerosol/dust in the air 
under the high dust. This implies a greater dust optical 
depth and therefore a greater chance for light scattered into 
the dust to be scattered back out toward the sensor. 
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Scattering of visible light by dust vice absorption is 
consistent with its low imaginary index of refraction at .86^m 
as indicated by Figure 2. The Temp 3 brightness temperatures 
are anomalously high and the reasons are in line with those 
discussed previously for "dust". 

c. Smoke: 

Smoke ranges from approximately 2.3 to 6.0 percent 
albedo on the Cal 2 axis and from 287°K to 306. 7°K on the temp 
3 axis. These ranges are as expected. The reasons for this 
signature characteristic were discussed in the Upper Gulf 
analys i s . 

d. Smoke/Dust Mixed: 

Smoke/dust mixed ranges from approximately 9.3 to 
20.0 percent albedo on the Cal 2 axis and from 312. 8°K to 
320®K on the Temp 3 axis. This signature exhibits the 
character istics of both smoke and dust. However, because of 
the relatively high albedos and temperatures, it appears the 
signature is dominated by the dust and/or the desert sand 
surface . 

e . Summary: 

This scatter plot emphasizes the significant 
difference in reflectance character istics between smoke and 
dust. Not only does smoke have a much lower albedo but its 
increase in channel 3 brightness temperature is not as 
linearly linked to increasing channel 3 reflectance component. 
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Dust, on the other hand, has a relatively high albedo and its 
high channel 3 brightness temperature is almost strictly 
influenced by the large channel 3 reflective component. The 
relationship is at least linear if not exponential. That is 
difficult to discern, however, because the AVHRR channel 3 
wavelength band upper limit restricts the detection of 
features warmer than about 322°K. That limitation accounts 
for the "bunching" of pixels vertically at 322°K on the 
scatter plot. Overall, the scatter plot distinguishes albedo 
signatures well but radically exaggerates "blackbody" 
temperatures. Additionally, even though clouds are not 
strictly analyzed, this image combination gives a clear 
separation between clouds and dust (high and low). 

2. Case 2 Scatter Plot Analysis: Cal 2 versus Temp 4 

Figure 25) 

There are various features in the Arabian Peninsula 
area (ie., dust, desert sand, clouds, etc.) that have similar 
solar reflectance characteristics. This image combination is 
used to discriminate between those features by its blackbody 
temperature . 

3. Dust: 

Dust ranges from 25.5 to 30.2 percent albedo on the 
Cal 2 axis and from 283®K to 288. 8®K on the Temp 4 axis. The 
Cal 2 range is in agreement with case 1. The Temp 4 range is 
significantly different from the Temp 3 range for dust (33®K 
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Figure 25. Case 2 Arabian Peninsula scatter plot. Feature 
designation corresponds with sample areas in Figure 23. 
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to 39°K difference). As indicated in Figure 2, dust has an 
order of magnitude higher imaginary index of refraction at 
10.8)im than at 3.7)im. That essentially makes it more of a 
"blackbody" in channel 4 than in channel 3. For this reason 
the Temp 4 range is more representative of the temperature at 
the top of the lower dust plume. 

High dust: 

High dust ranges from approximately 26.5 to 34.0 
percent albedo on the Cal 2 axis and from 277. 7°K to 283. 5°K 
on the Temp 4 axis. The Cal 2 range is consistent with case 
1. Once again, the Temp 4 range is indicative of the true 
temperature of the high dust plume as compared to the Temp 3 
range. The cooler temperatures compared to the lower dust are 
primarily altitude related. As previously stated, dust 
absorbs 10.8)im wavelength energy very well. Referring to 
Figure 26, it indicates the sun radiates little to no 
radiation at the 10.8)jlm wavelength. Therefore, the plume top 
dust is getting most of its 10.8)im energy from the atmosphere 
and so its temperature is representative of the atmospheric 
temperature at its particular altitude. Note that even though 
10.8)im is in a water vapor window, the sounding from Riyadh, 
Figure 6 shows a very dry atmosphere. This translates to 
little contribution by water vapor thermal emission to 



65 




Bmox 




Figure 26. Normalized Planck curves representing solar radiation 
( 5800K) and terrestrial radiation (260K) by wavelength (Jim), 
adapted from Kidder and Vender Haar, 1992) 
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EMITTANCE 



atmospheric temperature, further justifying Temp 4 plume top 
temperatures as indicative of atmospheric temperature. 

c. Smoke: 

Smoke ranges from approximately 2.3 to 6.0 percent 
albedo on the Cal 2 axis and from 283. 5®K to 290. 8®K on the 
Temp 4 axis. The Cal 2 range is consistent with case 1 and 
the Upper Gulf smoke albedo discussion applies. The Temp 4 
range is representative of thin smoke over land/water as 
indicated in the Upper Gulf cases. The highest temperatures 
represent thin smoke over water and the lower and temperatures 
represent thin smoke over land. The reasoning is also the 
same as in the Upper Gulf analyses. 

d . Smoke/Dust Mixed: 

Smoke/dust mixed ranges from 9.3 to 20.0 percent 
albedo on the Cal 2 axis and from 291. 5”K to 299®K on the Temp 
4 axis. The Cal 2 range is consistent with case 1. The Temp 
4 range is most representative of actual temperatures as 
compared to Temp 3. The signature discussion for albedo and 
temperature from case 1 is applicable here. 

e • Summary: 

Overall, this scatter plot depicts excellent 
signature separation between the four analyzed features. A 
final footnote, because Temp 4 has little reflective component 
in its brightness temperatures, the dust signature is almost 
all due to thermal emission and therefore most indicative of 
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its actual blackbody temperature and for reasons discussed 
previously, that of the atmosphere. There is, however, some 
ambiguity between the dust and cloud signatures as can be seen 
by the overlapping ovals on the scatter plot (Figure 25). 
This indicates an unambiguous representation of the dust 
distribution cannot be achieved when clouds are present using 
this image combination. 

3. Case 3 Scatter Plot Analysis: Cal 2 Versus Temp 5 

(Figure 27) 

This parameter combination was selected, primarily to 
investigate the difference in feature thermal emission between 
Chanels 4 and 5. All Cal 2 ranges are the same as previous 
cases . Therefore, only the Temp 5 signature ranges will be 
identified/ discussed. 

<3 . Dus t : 

Dust ranges from approximately 283. 8°K to 289. 8°K 
on the Temp 5 axis. These temperatures are slightly higher 
(~1°K) than in Temp 4 which alludes to an important dust 
characteristic. It implies dust is a better absorber /emi tter 
of channel 4 IR wavelengths (10.8pm). This is substantiated 
by Figure 3 from which a higher imaginary index of refraction 
can be inferred for channel 4 compared to channel 5. Because 
of this difference in absorption, the dust is not as optically 
thick in channel 5 compared to channel 4 and thus channel 5 
detects more thermal emmission from the warmer surface. This 
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accounts for the slightly higher channel 5 brightness 
temperatures. Additionally, note in Figure 28' that channel 4 
lies more in a water vapor window than channel 5. Therefore, 
channel 5 brightness temperatures are more sensitive to water 
vapor thermal emissions. Sounding data previously referenced, 
however, indicates a relatively dry atmosphere over the 
region; so the water vapor influence is most likely small. 
That is primarily speculation, however, because it is based on 
such sparse sounding data. 

b. High Dust: 

High dust ranges from 278. 2°K to 285°K on the Temp 
5 axis. Once again the temperatures are slightly higher 
( . 6°K to 1.5°K) than in Temp 4. The previous "dust" 
discussion as well as the discussion in case 2 for high dust 
apply here. 

c. Smoke: 

Smoke ranges from approximately 283. 5®K to 290. 8°K 
on the Temp 5 axis. This is consistent with the range 
analyzed in Temp 4, thus indicating no significant difference 
in smoke signature between the two channels. Therefore, the 
discussion in case 2 in applicable here. 
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Figure 27. Case 3 Arabian Peninsula scatter plot. Feature 
designation corresponds with sample areas in Figure 23. 
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Figure 28. Infrared water vapor (H 2 O) transmittance. Water 
vapor "window" is between 8-12jim. (Adapted from Kidder and 
Vender Haar, 1992.) The AVHRR/1 channel 4 and channel 5 
wavelength bands are delineated. 
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c?- Smoke/Dust Mixed: 



Smoke/dust mixed ranges from 291. 5°K to 299°K on 
the Temp 5 axis. This is the same range as in Temp 4 and that 
discussion from case 2 applies here. 
e . Summary: 

Overall, this scatter plot depicts the same 
information as in case 2. Of particular note again, is the 
ambiguity in signature between the dust and clouds. However, 
the one key difference is the higher temperature of the dust 
compared to Temp 4. This indicates dust is a slightly better 
"blackbody" in channel 4 than in channel 5. Therefore, 
channel 4 temperatures are probably more representative of the 
dust plume temperatures. 

4. Case 4 Scatter Plot Analysis: Cal 2 versus T4-T5 

(Figure 29) 

This parameter combination was selected to investigate 
the different dust absorption character istic between channels 
4 and 5. All Cal 2 ranges are the same as previous cases. 
Therefore, only the T4-T5 signature ranges will be identified/ 
d i scussed . 
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Figure 29. Case 4 Arabian Peninsula scatter plot. Feature 
designation corresponds to sample areas in Figure 23. 
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<3 • Dus t : 



Dust ranges from approximately -0.8°K to -1.47°K 
temperature difference on the T4-T5 axis. The negative 
difference for the dust signature is expected. The reasons 
for the difference were discussed in case 3. 

b. High Dust: 

High dust ranges from approximately -.35°K to - 
2.1°K temperature difference on the T4-T5 axis. The reasons 
for the negative differences were discussed in case 3. 

c. Smoke: 

Smoke ranges from approximately -3.2®K to +3.2°K 
temperature difference on the T4-T5 axis. The smoke is 
centered about zero which emphasizes its neutrality between 
channels 4 and 5. The positive difference is influenced by 
gulf water thermal emission and indicates channel 4's higher 
sensitivity to those emissions. The negative difference is 
influenced by the land thermal emissions and alludes to 
channel 5 being more sensitive to those emissions. 

d. Smoke/Dust Mixed: 

Smoke/dust ranges from approximately -1.8°K to 
+2.1°K temperature difference on the T4-T5 axis. This 
signature appeals to the neutrality of smoke and the albedo of 
a sand surface and/or dust. It is a complex signature and the 
location of the sample are on the scatter plot supports the 



hypothesis of a smoke/dust and sand surface signature 
combination . 

e . Summary: 

The primary significance of this scatter plot is 
that it emphasizes and quantifies the difference in dust 
signature between Temp 4 and Temp 5. The quantification 

allows for exploitation of the signature difference to locate 
dust outbreaks. Note on the scatter plot (Figure 29), 
however, there is a slight overlap between the dust and land 
areas 1 (Ll). This may pose some ambiguity in determining the 
dust distribution. 

5. Case 5 Scatter Plot Analysis: T3-T4 versus Cal 2 

(Figure 30) 

This image combination was selected to investigate the 
effect of the solar reflectance component in channel 3 on 
feature brightness temperature. The Cal 2 ranges are as 
analyzed in case 1. 
a. Dust: 

Dust ranges from approximately 32.5°K to 38.0°K 
temperature difference on the T3-T4 axis. The large positive 
difference values are indicative of the huge influence of the 

channel 3 reflective component on the brightness temperature 
of the dust. This large influence proves effective, however, 

in distinguishing the dust signature from the other features 
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Figure 30. 
des ignat ion 



Case 5 Arabian Peninsula scatter plot. Feature 
corresponds to sample areas in Figure 23. 
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when plotting it against its high albedo in channel 2. And 
unlike case 1, the upper channel 3 brightness temperature 
limit does not hinder the scatter plot from representing the 
true relationship between the reflective component of the 
channel 3 brightness temperatures and the Cal 2 albedos. 

b. High Dust: 

The high dust ranges from approximately 38.0°K to 
43.0°K temperature difference on the T3-T4 axis. As with the 
lower "dust", the large temperature difference influenced by 
the reflective component in the channel 3 brightness 
temperatures enhances the signature and creates excellent 
separation from the other features. 

c. Smoke: 

Smoke ranges from approximately 2.5°K to 16.0®K 
temperature difference on the T3-T4 axis. This is consistent 
with the ranges and discussion for thin smoke over land/water 
in the Upper Gulf analyses. 

d. Smoke/Dust Mixed: 

Smoke/dust mixed ranges from approximately 18.0°K 
to 25.0°K temperature difference on the T3-T4 axis. This once 
again is consistent with a signature representative of a 
combination of smoke and dust and/or sand surface. 

e - Summary: 

This scatter plot successfully separates the high 
dust and dust signatures from each other and the other 
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features. It is much more effective at showing the 
relationship between the channel 3 reflectance and channel 2 
albedo . 

6. Case 6 Scatter Plot Analysis: T4-T5 versus Temp 4 

(Figure 31) 

This image combination was selected because it gives 
excellent separation of the dust signature (both dust and high 
dust) from the other features on the scatter plot; the 
different absorption characteristic of dust between channels 
4 and 5 (as previously discussed) accounts for the negative 
dust values on the T4-T5 axis. There are some areas of land 
that have negative values but because those areas are warmer 
than the dust, they can be separated by plotting T4-T5 versus 
their temperatures (ie. Temp 4), as in this case. There are 
also some clouds that have negative to small positive T4-T5 
differences. Again, however, there is separation from dust 
based on temperature differences. Notice on the scatter plot. 
Figure 31, there is a low pixel density area (a line) running 
diagonally from left to right approximately from the point 
(275®K, +.4°K T4-T5 difference) to (303°K, -2.0°K T4-T5 
difference). All pixels below that area are dust and all 
other features are above the area. Thus, the conclusion is 
this scatter plot/image combination provides the most 
unambiguous representation of the dust signature. 
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Figure 31. Case 6 Arabian Peninsula scatter plot. Feature 
designation corresponds to sample areas in Figure 23. 
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c. VERIFICATION OF SCATTER PLOT ANALYSES 



As previously described in Chapter III, a two dimensional 
mask routine is utilized to verify the analyzed signature 
limits from the scatter plots. The mask was run on all cases 
but only the most significant cases from each area will be 
discussed . 

1. Upper Gulf Scatter Plots 

Two cases, S12L versus Temp 4 and T34 versus Temp 4, 
are presented for the Upper Gulf area. In both cases, only 
the smoke and water limits are masked. The thick smoke is in 
red, thin smoke over land in green, thin smoke over water in 
magenta and water in blue. 

a. S12L versus Temp 4 

Figure 32 is the scatter plot showing the boxed 
limits for the masked features. The limits correspond with 
those analyzed using the scatter plot technique. The limit 
selection is somewhat subjective, but notice that in this case 
there is relatively clear separation between features. S12L 
contrasts the smoke from the water and Temp 4 separates the 
smoke plume temperature from surface features and emphasizes 
the influence of surface emittance on the thin smoke radiative 
signatures (ie., as the smoke thins it begins to take on the 
signature of the surface over which it lies). Figure 33 is 
the mask of the limits for this scatter plot overlaid on the 
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Figure 32. Analyzed 2D feature radiative signature limits for 
case 1 Upper Gulf scatter plot. Colors correspond to those in 
mask. Figure 33. 
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S12L satellite image. Notice the thick smoke is primarily in 
the center of the plume and the thin smoke is on the periphery 
and downwind. The mask does a good job of representing the 
smoke distributions over both land and most importantly, 
water . 

b. T3-T4 versus Temp 4 

Figure 34 on the T3-T4 versus Temp 4 scatter plot 
shows the boxed limits for the masked features. (Note: thin 

smoke over land is described by two limit boxes.) The second 
limit was added to further delineate the feature. It is 
represented in range on the masked image. This image 
combination gives good separation of the main features - water 
and land from the smoke. On the other hand, it also 
emphasizes the complexity/ambiguity of the smoke signature due 
to the surface type influence. This is indicated by the 
overlap of the smoke limits. Figure 35 is the mask of the 
limits overlaid on the S12L satellite image. Comparison with 
the previous mask. Figure 33, shows the overall smoke 
distribution to be approximately the same. The significant 
difference is between the thick smoke versus thin smoke 
distributions. Those limits are more ambiguous on the T3-T4 
versus Temp 4 scatter plot and therefore are not as clearly 
masked as in the S12L versus Temp 4 case. Overall, however, 
it is a good representation of the smoke distribution. 
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Figure 33. S12L versus Temp 4 mask utilizing limits in Figure 
32. Mask overlaid on S12L image shows smoke distribution on 
1 March 1991 at 1026 UTC. Red is thick smoke, green is thin 
smoke over land, magenta is thin smoke over water and blue is 
water . 
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Figure 34. Analyzed 2D feature radiative signatures for Case 
3 Upper Gulf scatter plot. Color corresponds to those in 
mask. Figure 35. 



84 








Figure 35. T3-T4 versus Temp 4 mask utilizing limitations ii 
Figure 34. Mask overlaid on Temp 4 image shows smok< 
distribution on 1 March 1991 at 1026 UTC . Red is thick smoke 
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2 . 



Arabian Peninsula Scatter Plots 



Three cases are discussed for this area - Cal 2 versus 
Temp 5, Cal 2 versus T4-T5 and T4-T5 versus Temp 4. Only the 
dust limits are masked for each case. 

a. Cal 2 versus Temp 5 

Figure 36 is the scatter plot depicting the limit 
boxes for low dust and high dust. The difference between low 
dust and high dust is based on temperature with low dust being 
warmer. The mask of these limits. Figure 37, represents low 
dust in red and high dust in yellow with the Temp 5 satellite 
image underlaid. The mask, however, is somewhat ambiguous in 
representing the dust areas. Referring to the Cal 2 (Figure 
38) and Temp 5 (Figure 39) satellite images, some of the 
clouds, as well as some of the coastal land areas along the 
Red Sea and Arabian Sea, have the same signature as the dust 
in both images; they are masked along with the dust. 
Therefore, even though the dust plume is represented, it is 
not done so unambiguously. 

b. Cal 2 versus T4-T5 

Figure 40 is the scatter plot depicting the limit 
boxes. In this case, the dust limits were subdivided into 
four regions. This gives better resolution and depiction of 
the dust distributions. Referring to the Cal 2 (Figure 41) 
and T4-T5 (Figure 42) satellite images, it can be seen parts 
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Figure 36. Analyzed 2D ”dust” radiative signature for Case 3 
Arabian Peninsula scater plot. Colors correspond to those in 
mask. Figure 37. 
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Figure 38. Cal 2 1 March 1991 1026 UTC 
image. Some clouds (circles) and coastal 
have similar brightness as dust plume. 



Arabian Peninsula 
land (rectangles) 



89 




Figure 39. Temp 5 1 March 1991 1026 UTC Arabian Peninsula 
image. Some clouds (circles and coastal land (rectangles) 
have similar brightness temperatures as dust plume. 
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Figure 40. Analyzed 2D "dust" radiative signature for case 4 
Arabian Peninsula scatter plot. Colors correspond to those in 
mask. Figure 41. 
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of the boxed land area has the same signature as the dust in 
both images. As in the previous case, this causes some 
ambiguity with the mask. Figure 43 shows the mask overlaid on 
the Cal 2 satellite image. The lowest (warmest) dust is in 
red followed by yellow, green and the highest (coldest) dust 
in blue. Notice, in fact, much of the boxed land area is 
masked along with the dust. Thus, once again the dust 
distribution is depicted but not clearly. 
c. T4-T5 versus Temp 4 

Figure 44 shows the boxed limits for the dust. 
Because the dust was clearly separated from the other features 
it was subdivided into smaller limits so the mask can give 
better resolution on its distribution. Figure 45 is the mask 
of this image combination overlaid on the Temp 4 satellite 
image. Low (warmest) dust is in red followed by orange, 
yellow, green, blue and the highest (coldest) dust in magenta. 
The highly discernible dust limits allow this mask to be the 
most unambiguous representation of the dust plume. 
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Figure 41. Cal 2 
Boxed land area has 



1 March 
s imi lar 



1991 1026 UTC 

albedo as parts 



satellite image, 
of dust plume. 
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Figure 42 . 
Boxed land 



T4-T5 
ara has 



1 March 1991 1026 UTC satellite image, 

similar T4-T5 difference as dust plume. 
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Figure 43 
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CH4 TEMP. (K) (10.8-mM) 



Figure 44 . Analyzed 2D "dust" radiative signature for Case 6 
Arabian Peninsula scatter plot. Colors correspond to those in 
mask. Figure 45. 
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p£gujr 0 45. T4""T5 VS 3 TSUS T6inp 4 inssk utilizin <3 lirnits in 
Figure 44. Mask overlaid on Temp 4 image shows unambiguous 
dust distribution. Where red represents the lowest dust and 
magenta represents the highest dust. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. SMOKE RADIATIVE SIGNATURE 

The following conclusions have been drawn concerning the 
smoke radiative signature of smoke. First, smoke detection 
over water is more difficult than over land. Isolation of the 
signature is complicated by the low reflectance of both smoke 
and water in the solar spectrum and relative transparency of 
smoke at IR wavelengths. The image combination of S12L versus 
Temp 4 is the best for discerning smoke over water. 
Contrasting the different absorption characteristics of water 
between channels 1 and 2 isolates the smoke signature. 
Second, smoke absorbs the visible spectrum wavelengths equally 
well (appears black on all visible images). Thus, there is 
strong contrast with most land regions which makes it 
relatively easy to detect. Third, smoke particles become more 
transparent with increasing IR wavelength. This is apparent 
from image combination of T3-T4 versus Temp 4. Fourth, high 
atmospheric water content in the presence of smoke particles 
provides better IR wavelength detection. The smoke particles 
act as condensation nuclei and allow liquid water drops to 
form. And finally as smoke thins, its signature becomes more 
and more characteristic of the underlying surface which 
complicates its detection as it disperses. 
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B. DUST RADIATIVE SIGNATURE 



There are four main conclusions drawn in this study 
concerning the dust signature. First dust detection over land 
is more difficult than over water. Isolation of the aerosol 
signature is complicated by the influence of the highly 
diverse land features. Second, dust is highly reflective at 
visible wavelengths. Third, the dust exhibits a general trend 
of increasing absorption with IR wavelength. Finally, the 
strong difference of absorption between channels 4 and 5 
allows for the most accurate detection of dust. 

C. TWO DIMENSIONAL MASK AND SCATTER PLOT TECHNIQUE 

The mask routine uniquely illustrates aerosol spatial 
distributions using two dimensional radiative limits. It 
clearly defines detected aerosols in a way that does not "wash 
out" when reproduced or copied. Since it represents aerosol 
distributions well, it would be an excellent mission briefing 
t ool . 

The scatter plot technique works well in separating and 
defining aerosol radiative signature limits. Of course, some 
image/parameter combinations are better than others. But the 
idea is to find the best combination that contrasts with known 
characteristics of the aerosol and separates it from the other 
features as was done in this project. 

The scatter plot technique in combination with the mask 
routine does show promise for possible fleet application. 
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Granted, there is much more work to be done, but it is 
conceivable that an algorithm based on these techniques could 
one day be used for aerosol tracking and mission 
planning /briefing. 

D . RECOMMENDATI ONS 

The two parameter scatter plot and mask, in certain cases, 
gave some relatively unambiguous depictions of aerosol 
distributions. Investigation of a three parameter scatter 
plot/mask routine, however, is recommended. Three parameters 
should further enhance aerosol signature detection. As an 
example, unique smoke characteristics are represented in S12L, 
Temp 4 and T3-T4 parameters. Analysis in all three, 
simultaneously, would allow even more unambiguous signature 
separation than is currently attained relying on just two 
parameters . 

It is recommended to build upon the combination of the 
scatter plot analysis technique and masking routine to devise 
an aerosol detection algorithm. As implied previously, such 
an algorithm could have military application. 

This project was conducted using a single day case study. 
The masking limits analyzed for the aerosols in this project 
are most likely only valid for 1 March 1991. Therefore, it is 
recommended that additional cases be examined to further 
verify the techniques used in this project and to study 
geometry and temperature variation effects. 
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Finally, it is recommended that future AVHRR radiometers 
have a higher upper end channel 3 wavelength limit. The 
current limit does not allow "hot" objects greater than 
approximately 322“ K to be discerned. This is a result of not 
accounting for the effect of solar reflectance on channel 3 
brightness temperatures. 
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